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Abstract. In this paper we analyse and discuss the HR Diagram position of two extra-solar planet-host stars - 
HD 37124 and HD 46375 - by means of theoretical stellar evolution models. This work was triggered by the results 
obtained by Laws et al. 1120031 who found that these stars were in contradiction to the expectation based on their 
high metallicity. 

Fixing the age of both stars with the value based on their chromospheric activity levels and computing our 
own evolutionary models using the CESAM code, we are able to reproduce the observed luminosity, effective 
temperature and metallicity of both stars for a set of stellar parameters that are astrophysically reliable even if 
it is non-trivial to interpret the absolute values for these parameters. Our results are discussed in the context of 
the stellar properties of low mass stars. 
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1. Introduction 

During the last few years a number of different de- 
tailed spectroscopic studies have revealed that planet- 
host stars are (on average) very metal-rich compared to 
field dwarfs (e.g. Gonzalez et al. 1953 IWlTl Santos et 
al. I^nroi IMTTl HUM Laws et al. gOOSjl . Current results 
favour that this "excess" metallicity arises from the cloud 
of gas and dust that generated the star/planetary-system, 
suggesting that the formation of giant planets (at least 
the kind currently beeing discovered) is critically depen- 
dent on the grain content of the proto-planetary disk (e.g. 
Santos et al. I2(JU4|I . This result has an enormous impact 
on theories of planetary formation. 

In one of the studies mentioned above, Laws et 
al. H2U(J3 1 analysed a sample of 30 F, G, K stars 
with companions (planets or brown dwarfs). They pre- 
formed spectroscopic observations in order to determi- 
nate their atmospheric parameters, including the effec- 
tive temperature and metallicity. All the stars had accu- 
rate fflPPj4iZC0S' parallaxes, and thus the stellar position 
in the Hertzsprung-Russel Diagram (HRD) could be ob- 



tained precisely enough to derive stellar masses and ages 
from theoretical stellar models. To derive these parame- 
ters. Laws et al. (2003) used the evolutionary tracks from 
the Padova group (Salasnich et al. 2000 ) . In their conclu- 
sions it is pointed out that for four stars in their sample 
(HD6434, HD 37124, HD 46375, HD 168746) the derived 
ages were above 14.8 Gy in spite of being Population I 
stars. They called this quartet the "red stragglers", be- 
cause of their unexpected red position in the HRD. 

The explanation for the position of these four stars in 
the HRD is not clear. Laws et al. f20g3|l suggested that 
it could be related to the presence of close-in low-mass 
stellar companions. However, to our knowledge none has 
been found by any Adaptative Optics survey. Laws et al. 
H2003|l "recommend that these stars be the targets of de- 
tailed studies to determine the source(s) of their aberrant 
physical characteristics" . We perform this study starting 
from a detailed analysis of the position of two of these 
objects in the HRD, to verify to what extent these two 
objects can be explained by current models of stellar evo- 
lution. 
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The knowledge of the HRD position of a star is depen- 
dent on many different variables. To first order, it is fixed 
by the stellar mass (A/), the initial individual abundance 
of helium (Y) and metals {Z), and age {t^,). The resolu- 
tion of the equations of internal structure using the above 
parameters will give the values for the surface tempera- 
ture and luminosity. However, the physical inputs chosen 
to describe the stellar interior also constrain the evolution 
of the star in the HRD. In particular, some mechanisms 
insufficiently known such as the convection, rotation and 
diffusion are dependent on free parameters. For instance, 
in the framework of the Mixing Length Theory (MLT, 
see below) currently used to model stellar convection, one 
more (unknown) variable has to be considered: the mixing- 
length parameter (a). So, to model a star correctly, we 
must determine the five parameters. 

From the observational point of view, a detailed spec- 
troscopic analysis of a single star may allows us to estimate 
its metalhcity {[Fe/H] or Z/X, where X = 1~Y-Zis 
the abundance of hydrogen), its luminosity (L)^ and its 
effective temperature. Teg. We are thus left with three 
known observables. In other words, we have an indetermi- 
nate problem: five unknowns to three observables. Gravity 
is also an issue in the spectroscopic analysis. So, a deter- 
mination of the mass should, in principle, also be possible. 
However, a typical realistic error in logg is usually above 
0.1 dex {e.g. Laws et al. 2003). This error implies an un- 
certainty on the derived solar mass stars of more than 0.15 
Mq . Clearly, this error is too large to allow the gravity to 
be a real mass constraint. 

This indeterminate problem can be solved for the Sun, 
for which the mass, metalhcity and age are observationally 
known (Christensen-Dalsgaard,1982J and for some binary 
systems where individual masses are available {e.g. Noels 
et al. 119911 Fernandes et al. I1998|l . But for a single star 
(without a close-in stellar companion), other than the Sun, 
commonly there are more unknowns than observables, the 
situation that we have here. 

In this work we analyze in detail two of the "red- 
stragglers" discussed by Laws et al. I)2()()3|l . namely 
HD 37124 and HD 46375. Our goal is to find physically 
consistent models that are able to explain the position of 
these stars in the HRD. We have excluded HD 168746 and 
HD 6434 from this study: the former because it has no 
available chromospheric activity age (which would give us 
a further constraint - see below), and the latter because 
with a [Fe/H] = —0.55 ±0.07 it is on the frontier between 
Population I and Population II stars. In the modelling of 
Population II stars some effects become important. This is 
the case of the diffusion of chemical elements, NLTE cor- 
rections for the metalhcity and a-elements contribution to 
the opacity (Lebreton et al. I1999|l . So, as the cumulated 
effects can be a supplementary source of errors, we prefer 
to exclude this case. 



^ Better determined if the distance to the star is known, as 
it is in our case. 



In Sec. 2 we present the stellar models and the theo- 
retical description of our methodology used for modelling 
the stars. In Sec. 3 we compute models for HD 37124 and 
HD 46375. In Sec. 4 we present a discussion and conclu- 
sions in the context of the stellar proprieties of low mass 
stars. In particular, we show that the two stars studied 
in the current paper can be correctly modelled for some 
physically reliable values of the stellar parameters. 

2. Input physics in the theoretical stellar models 
and modelling method 

We compute our models specially for the stars studied 
in this work. The stellar evolution calculations were done 
with the CESAM code version 3 (Morel [TM7|l . running in 
the Coimbra Observatory, in Portugal. 

Details on the physics of these models can be found 
in Lebreton et al. 1^1999^1 : the CEFF equation of state, in- 
cluding Coulomb corrections to the pressure (Eggleton et 
al. I1973L Christensen-Dalsgaard ITI)91|I : nuclear reactions 
rates given by Caughlan and Fowler H1988|l : solar mixture 
from Grevesse and Noels (|1993|l ; OPAL opacities (Iglesias 
and Rogers [T996|l complemented at low temperatures by 
opacity data from Alexander and Ferguson fl993') follow- 
ing a prescription of Houdek and Rogl (^996). The atmo- 
sphere is described with an Eddington T(r)-law; convec- 
tion is treated according to the mixing-length theory from 
Bohm-Vitense p958f) leaving the mixing-length {a x Hp) 
unknown, and thus a as a free parameter, where Hp is the 
pressure scale height. With these input physics the solar 
model fits the observed luminosity and radius, with a rel- 
ative precision of 10"'^, with a=1.63, helium abundance 
y=0.268 and Z = 0.0175 for the common accepted solar 
age of 4.6 Gyr (Dziembowski et al. 11999(1 and the ratio of 
the solar mixture from Grevesse and Noels (1993). This 
model does not take into account microscopic diffusion of 
helium or metals. 

There are other mechanisms, not included in this work, 
that could affect the HRD position of a model, such as the 
rotation (Maeder and Meynet ,2p00) , overshooting of the 
convective core, or helium and metal gravitational set- 
tling (microscopic diffusion). However, the stars is this 
work are slow rotators (as are all planet-hosts, since it is 
more difficult to obtain accurate radial-velocities for fast 
rotating stars) and solar mass stars do not develop per- 
manent convective cores (Ribas et al. 1200011 . On the other 
hand, diffusion has a marginal effect on the HRD position 
for Population I solar mass stars. Lastennet et al. I|2()()3|l 
have shown that for a O.98M0 star with Z=0.012, the in- 
clusion of helium and metal diffusion may change the Tes 
of the evolutionary track by no more than 65 K, even for 
an old age of 10 Gy. For Population I stars, the density in 
external layers is sufficiently high to break the diffusion ef- 
ficiency. We thus expect that the lack of these mechanisms 
do not considerably change our results. 

As discussed above, modelling a single star, for a fixed 
age, is reduced to the problem of finding the solutions of 
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the four stellar parameters Af , Y, Z and a, that adjust 
the three observational parameters L, Toff, and Z/X . 
Inspired by solar modelling (Christensen-Dalsgaard 

EMI, 

we construct the following system composed of the 
three equations: 

L L dlog-r- 

(log— )o6. - i\0g—)ref+J2i^^)j^rXiX,^X,^Ji 

logTeff^,^ = logTeff„^+^(^i^i^),^.x(X.-X.„,)(2) 



where Xi with i = 1,2,3, and 4 are Af, y, Z, a, respec- 
tively. The subscript "ref" indicates the reference stellar 
evolutionary model that falls inside the observational er- 
ror bars. We consider as a satisfactory solution the sam- 
ple (Af, Y,Z, a), obtained by solving the above system 
assuming the linearity of the variation of the logL/LQ 
and log Tcff with mass, helium, metals and a for a fixed 
value of the age. Theoretically, each solution of the system 
should fit exactly the observed luminosity and effective 
temperature, for the fixed age. In practice, thanks to our 
linear approximation, in some rare situations the solution 
is outside the error bars. However, in these situations the 
solution is still near and a small change in one parame- 
ter is enough to put the evolutionary track in the error 
bar. On the other hand, as we have more unknowns than 
equations, the above system has an infinite number of so- 
lutions. In this work we chose to present five solutions for 
each star in order to illustrate the range of variation of 
the stellar parameters (Af, F, Z, a). 

3. Modelling HD 37124 and HD 46375 

In the following we will use the ages for HD 37124 
and HD 46375 that have been estimated by Laws et al. 
I|2()()3|l using a chromospheric activity-age calibration from 
Donahue 1)1993(1 . The knowledge of this age gives us the 
possibility to limit the number of models available. We 
have thus fixed the age for each star (3.9 and for 4.5 Gyr 
for HD 37124 and HD 46375, respectively) and estimated 
the other four unknown parameters (A/, y, Z, a) by theo- 
retical stellar evolutionary models, taking into account the 
three observables (L, Teff ,Z/X). So, opposite to Laws et 
al. (2003), we have chosen not to fix the values of a and 
in our models. We note that our major goal is to 
find an astrophysically reliable solution for the position 
of these stars in the HRD, i.e. a solution that matches 
the observable quantities. We thus prefer to be as strict 
as possible in our observable constraints. Leaving the age 
as a free parameter would, of course, make this solution 
easier to obtain, but this solution would probably not be 
reliable. 



Table 1. Observations of HD 37124 and HD 46375 (see 
text for details). 



star 


L/Le 


Teff (K) 


Z/X 


HD 37124 
HD 46375 


0.84 ±0.06 
0.74 ±0.06 


5551 ± 34 
5241 ± 44 


0.0105 ± 0.0008 
0.049 ± 0.004 



3.1. HD 37124 

3.1.1. HD 37124: observations 

HD 37124 is a G4 IV-V star (V==7.68; B-V==0.667) 
in the solar neighbourhood at a distance of 33.2 pc 
{HIPPARCOS -Peiryman et al. lTWTjl . Vogt et al. liM17)|l 
reported the discovery of a candidate planet around this 
star, followed by the discovery of a second planetary-mass 
companion in the system (Butler et al. l2003;l . Since then, 
particular attention has been paid to this object, e.g., de- 
tailed spectroscopic analyses (Laws et al. I20U3I Santos et 
al. I2004|l or studies of the dynamical stability of the dis- 
covered planetary system (Zhou and Sun l2003p . 

In the following we will consider the spectroscopic de- 
terminations and magnitudes given by Laws et al. plH^ : 
Toff =5551 ± UK; [Fe/H] = -0.37 ±0.03 and My = 
5.07 ± 0.08. These spectroscopic results are in excellent 
agreement with the recent values obtained by Santos et 
al. fSOOl: Toff =5546 ±30i^; [Fe/H] = -0.38 ± 0.04. 

For Population I stars, the abundance of metals Z is re- 
lated to the [Fe/H] by [Fe/ H]=log{Z/ X)-\og{Z/ X)q, 
where (Z/X)q=0.0245 is the ratio of the solar 
mixture (Grevesse and Noels I1993|l . This implies 
(Z/X)=0.0105±0.0008 for HD 37124. 

Using the T^g of Laws et al. 1)2003(1 we derive a bolo- 
metric correction from Flower 1)1996(1 of BC — —0.13 ± 
0.01. This implies a bolometric magnitude of Mboi = 
4.94 ± 0.08. Assuming that the bolometric magnitude of 
the Sun is 4.75 (e.g. Gayrel I2002p . we obtain a lumi- 
nosity of 0.84±0.06 times the solar luminosity {Lq) for 
HD 37124. From the Stefan-Boltzmann law we can thus 
estimate the radius of this star: 0.99±0.05 Rq. This result 
is in very good agreement with the infrared photometry 
determination of 1.004 ± 0.039 from Ribas et al. l(M7Sjl . 
The observable stellar parameters used to constraint the 
models are summarised in Tabled 

3.1.2. HD 37124: models 

We have computed several models for HD 37124, fixing 
the age at 3.9 Gy from Laws et al. (20®, in the follow- 
ing range of parameters: 4 values of mass in [0.75, 0.95], 4 
values of helium in [0.24,0.30], 4 values of metallicity in 
[0.00669, 0.01000] and 5 values of a in [0.8, 1.6]. This sam- 
ple of models allows us to obtain the derivatives of L and 
Tcff in relation to (Af, Y, Z, a) in the system of equations 
(1), (2) and (3). Using the observed values (T, Tcs, Z/X) 
we can then establish a linear system for the reference 
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Table 2. Models of HR 37124 for an age of 3.9 Gy. M 

is the mass of the model (in Mq), Y is the initial helium 
abundance, Z is the initial metal abundance, a the mixing 
length (in units of the local pressure scale height). 



Model 


M 


Y 


Z 


a 


Toff (K) L/Lq 


model 1 


0.86 


0.29 


0.0074 


0.85 


5548 


0.84 


model 2 


0.90 


0.27 


0.0076 


0.90 


5526 


0.84 


model 3 


0.92 


0.25 


0.0079 


0.95 


5535 


0.84 


model 4 


0.94 


0.24 


0.0079 


1.00 


5542 


0.84 


model 5 


0.90 


0.27 


0.0076 


1.10 


5668 


0.84 



why Laws et al. (2003) did not find a satisfactory solution 
in their analysis^. In fact, all the Padova isochrones are 
built using the solar a value 1.7), the same for all stel- 
lar models. This approximation, currently made on the 
isochrone computations for simplicity, is not clear: why 
should stars with different masses, chemical composition 
and evolutionary status, have the same a ? Some recent re- 
sults on stellar systems conclude that there are variations 
of a with mass (Lebreton et al. I2U(J2|I . including values 
of a lower than 1.0 (Lastennet et al. l2U(j3|l . On the other 
hand, the helium and mass are values expected for this 
star. 



.J 



-0.05 



-0.1 



-0.15 



-0.2 




HD37124 



3.75 



J I I I I , I I 

3.745 3.74 
log(T,„) [k] 



J I L 

3.735 



Fig. 1. Evolutionary tracks for HD 37124. The different 
lines represent models 1 (solid line), 2 (dotted line), 3 
(dashed line), and 4 (dot-dashed line). See text for more 
details. 



model with the following characteristics: M=0.94 Mq, 
F=0.24, 2'=0.0079, and a=1.00. This model reproduces 
the observed quantities (L, Tcs , Z/X), within the errors, 
for an age of 3.9 Gyr. 

In Table [3 we present the models from 1 to 4 that 
reproduce the observations within the error bars and in 
Fig. ^ we provide the corresponding evolutionary tracks 
in the HRD from the ZAMS to 3.9 Gy. The number of 
possible models that fall inside the observed error bars 
is infinite. The first four models were chosen in order to 
show the dependence of the solution on M, Y, Z and a. 
Model 5 is computed taking into account an increase in the 
observed Tes of about 3a (see Discussion and conclusions). 
The correponding luminosity is 0.82^0, still inside the 
original error bar. 

The five solutions have in common low values of the 
mixing-length parameter, a. This is probably the reason 



3.2. HD 46375 

3.2.1. HD 46375: observations 

HD 46475 is a Kl IV-V star (V=7.91; B-V=0.860), lying 
in the solar neighbourhood: d = 33.4 pc fPerrvman ll997fl . 
This star was recently found to be orbited by a possible 
planetary-mass companion (Marcy et al. I2()()()|l . 

To build the observed HRD we consider here the pa- 
rameters derived by Laws et al. (2003): Tcff =5241 ±44/^; 
[Fe/H] = +0.30 ± 0.03 and My = 5.29 ± 0.08. These 
spectroscopic results are in good agreement with the re- 
cent results by Santos et al. (2004) for the effective tem- 
perature (Toff =5268 ± 55 K) and are in fair agreement 
with the metallicity {[Fe/H] = -t-0.20 ± 0.06). Computed 
as described above, the metal abundance of this star is 
(Z/X)=0.049±0.004. 

Using the Toff of Laws et al. H2003fl we can derive the 
bolometric correction of = -0.21±0.03 fFlower lTU^ . 
a value that implies a bolometric magnitude Mboi = 5.08± 
0.09, and a luminosity of 0.74±0.06 L©; the corresponding 
stellar radius is 1.05 ± 0.06 Rq, once more in very good 
agreement with Ribas et al. pUUHIl : 1.005 ± O.O36i?0. The 
observations chosen to constrain the models are provided 
in Tabled 



3.2.2. HD 46375: models 

As done for HD 37124, we computed several models for 
HD 46375, fixing the age at 4.5 Gyr from Laws et al. l^dU^ 
in the following range of stellar parameters: 4 values of 
mass in [0.95, 1.05], 4 values of helium in [0.24, 0.30], 4 val- 
ues of metals in [0.025, 0.035] and 4 values of a in [0.8, 1.6]. 

The reference model has the following chacteristics: 
Af=1.00 Mq, r=0.30, Z=0.036 and a=1.0. This model 
is sufficiently near the observed quantities {L, T^ff , Z/X) 
for an age of 4.5 Gyr. 

In Table |31 we present four derived models, from 6 to 9, 
that reproduce the observations within the error bars and 
in Fig.|5]we plot the correspondent evolutionary tracks in 
the HRD from the ZAMS to 4.5 Gyr. As for the previous 
star we chose four models to give the range of variation 



^ By "satisfactory solution" we mean the one that, for a fixed 
age, reproduces the observational constraints (L, T^s, Z/X). 
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Table 3. Models of HD 46375, for an age of 4.5 Gy. M 
is the mass of the model (in Mq), Y is the initial helium 
abundance, Z is the initial metal abundance, a the mixing 
length (in units of the local pressure scale height). 



Model 


M 


Y 


Z 


Q 


Teff {K) L/Lq 


model 6 


0.96 


0.31 


0.032 


1.15 


5236 


0.71 


model 7 


0.98 


0.31 


0.032 


1.05 


5224 


0.79 


model 8 


1.01 


0.29 


0.033 


1.10 


5200 


0.77 


model 9 


1.03 


0.27 


0.034 


1.35 


5272 


0.71 


model 10 


1.03 


0.27 


0.034 


1.55 


5369 


0.72 



-0.1 



-0.15 



.J 

M no 
o —0.2 



-0.25 



-0.3 




HD46375 



3.725 



3.72 3.715 
log(Tj [k] 



3.71 



Fig. 2. Evolutionary tracks for HD 46375. The different 
lines represent models 6 (solid line), 7 (dotted line), 8 
(dashed line), and 9 (dot-dashed line). See text for more 
details. 



of each parameter Af , Y, Z and a. Model 10 is computed 
taking into account an increase of 3cr in the observed Tee 
. As previously, the corresponding luminosity is O.7IL0, 
falling inside the original error bars. 

As for HD 37124, the models from 6 to 9 present low 
a values confirming that a probable explanation for the 
difficulty of Laws et al. H2U(J3|I in finding a solution with a 
"reliable" age can be attributed to the constant value of 
a in the Padova models. 



4. Discussion and conclusions 

In this paper we present detailed evolutionary models that 
reproduce the observed HRD position and metallicity of 
the planet-host stars HD 37124 and HD 46375. In order to 
reduce the high degeneracy degree of modelling a single 



star we fixed the stellar age given by the age-cromospheric 
activity relation. 

The main goal of this paper is to show that it is pos- 
sible to model these stars using astrophysically reliable 
values of mass, helium and metal abundances, and mix- 
ing length parameters. We found that this is indeed the 
case, although a and values are not considered as 
a constants from one model to another. The hypothesis 
of constant values is current in published isochrones. We 
propose that the "aberrant" solutions found by Laws et 
al. [ 20W) for HD 37124 and HD 46375 are caused by the 
assumption of the above hypothesis, particularly the one 
about the mixing length parameter. 

Eventual differences between the internal structure of 
the CESAM and Padova evolutionary codes, as well as 
between our modelling method and the one used by Laws 
et al. H2Ul)3|l . can't explain our results. The three follow- 
ing reasons support this: first, the input physics, namely 
the radiative opacities, equation of state and convection 
formalism, are similar in the two codes; second, we con- 
firm the difficulty of Laws et al. I|2()()3|l in modeling the 
above stars using solar values for a and ; and third, we 
tested our modelling method in HD 117176, one star also 
studied by Laws et al. I|2U(J3|I . and for which they have 
found a satisfactory solution. Using the mass derived by 
these authors, as well as the same a and used in the 
Padova isochrones, we find a solution in the HRD for an 
age of 8.8 Gy, close to the 7.9 Gyr derived by them. 

Among an infinite number of combinations of stellar 
parameters (M, Y, Z, a) that, for a fixed age, reproduce 
the observables (L, Tcs , Z/X)^ we have computed four 
models for each star, inside an acceptable range of varia- 
tion for each parameter. As we do not have complemen- 
tary observational constraints we cannot choose the best 
model(s). However, some astrophysical arguments, such 
as the reliable values of the helium-to-metals chemical en- 
richment parameter (e.g. Jimenez et al. ''2003) could 
allow us to constrain the solutions. Note that our solutions 
are age dependent. 

Understanding why these stars are well modelled us- 
ing a low a value, while that is not the case for other 
stars in the sample of Laws, is beyond the scope of this 
paper. The MLT is a local, ID and non-compressive the- 
ory to explain a non-local, 3D and highly compressible 
system, as is stellar convection. Many works and authors 
(e.g. Canuto and Mazzitelli 1991 Ludwig et al. ll999|l . have 
already shown the limitations of this formalism. However 
MLT is easy to implement in stellar evolutionary codes 
and is a relatively successful in reproducing global obser- 
vations of the Sun and stellar systems. The low values of 
a found in this work, probably, reflect the limitations on 
the MLT itself more than a real physical effect. The anal- 
ysis of these stars using alternative convective description 
such as the formalism of Canuto and Mazzitelli II1991II 



As mentioned above, the modelling of a single star other 
than the Sun is an indeterminate problem, given that there are 
more unknowns than observables. 
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could give important constraints on the super-adiabatic 
layer in the convective region, where convective efRciency 
is strongly dependent on a. However, we cannot exclude 
that the derived low values for a are due to some other 
unknown effect. 

It is well known that the T^s of low mass stellar mod- 
els depends on a, thanks to its influence in the super- 
adiabatic layer and, thus, on stellar radius. Are there er- 
roneous Toff determinations? We have checked that if the 
observed Toff is increased by 3ct for each star, we can find 
a solution with a high a value for both stars (models 5 
and 10). In particular, model 10 has an a value closer to 
solar. On the other hand these changes in Toff must be 
weighted by the fact that systematic differences between 
various observational temperature scales are never above 
100-200 K (c.f. Laws et al. 12003 Santos et al. l^UMjl . So at 
least for HD 46375 the reason for the "abberant solution" 
found by Laws et al. (j2003|l may be an underestimation of 
Tcfi . 

The solutions found for each star are all in a range 
lower than 0.1 solar masses. This can be particularly im- 
portant for the estimation of the mass the planets har- 
boring stars, and that could have a strong impact e.g. on 
the study of the relation between the stellar mass and the 
frequency of giant planets (Santos et al. 120031 Laws et al. 
[2003J. Furthermore, the wrong masses derived from the 
standard evolutionary models (with, among other things, 
a fixed to the solar value) may imply incorrect derivations 
of the surface gravities using Hipparcos parallaxes (e.g. 
Nissen et al. I1997|) . If the values of a depend e.g. on the 
effective temperature of the star, this problem may induce 
a systematic error in the gravities with temperature when 
comparing spectroscopic and parallax-based \ogg values, 
an effect that was mentioned in Santos et al. H2U04|I . 

If indeed parameters like a vary from star to star, the 
determination of accurate ages by means of isochrone fit- 
ting may introduce important systematic errors. This was 
the case for the two stars studied in this paper, but may 
be the case for many other objects thus, specific stellar 
models must be computed for individual objects. 
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